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We introduce a single-channel opening in a random Hamiltonian and a quantized chaotic map: 
localization on the opening occurs as a sensible deviation of the wavefunction statistics from the 
predictions of random matrix theory, even in the semiclassical limit. Increasing the coupling to the 
open channel in the quantum model, we observe a similar picture to resonance trapping, made of 
few fast-decaying states, whose left (right) eigenfunctions are entirely localized on the (preimage of 
the) opening, and plentiful long-lived states, whose probability density is instead suppressed at the 
opening. For the latter we derive and test a linear relation between the wavefunction intensities and 
the decay rates, similar to Breit-Wigner law. We then analyze the statistics of the eigenfunctions 
of the corresponding (discretized) classical propagator, finding a similar behavior to the quantum 
system only in the weak-coupling regime. 

PACS numbers: 03.65.Ta,05.70.Ln,89.70.Cf,05.70.-a 


I. INTRODUCTION 

One of the distinctive traits of all chaotic systems is 
their seemingly ‘random’ behavior [l| . As a consequence, 
one usually assumes that the eigenfunctions of a quan¬ 
tized chaotic Hamiltonian have the same statistical prop¬ 
erties (j.e. wavefunction intensity distribution) of a com¬ 
plete set of waves with random amplitudes and phases 
13,3 , or equivalently, of the eigenvectors of a Hermitian 
matrix with random entries, according to random matrix 
theory (RMT) 0,0. Due to a number of applications 
(quantum information theory 013| , classical ja ^3 and 
quantum optics EMI, quantum transport [3,^0), as 
well as to equally many theoretical issues (see for exam¬ 
ple [313), quantum chaos community is nowadays 
largely focused on the behavior of open systems [l3|. 

In this paper we address one of the simplest theoreti¬ 
cal questions: whether and how the wavefunction statis¬ 
tics deviates from the predictions of the random wave 
assumption as we perturb a chaotic system with a single¬ 
channel opening. As main result of our investigation, 
we numerically find that the overall wavefunction in¬ 
tensity distribution at the location of the opening does 
change from the RMT-expected shape to a longer- 
tailed curve, which is analytically described using per¬ 
turbation theory. It physically implies that localization 
occurs at the opening. In our theory the opening can 
be an arbitrary state |a) in the Hilbert space, however, 
in most of our testing models we take it as a coherent 
state in the phase space. Deviations of the wavefunction 
statistics from RMT have been observed before in real 
space: for time- reversally symmetric systems, it was con¬ 
jectured [3| and then shown analytically and experimen¬ 
tally [2l| that the distribution of the wavefunctions at the 
leads smoothly crosses over from Porter-Thomas’ to Pois¬ 
son’s with the coupling to the opening. Although there 
was no explicit mention of localization, the wavefunction 


distribution for a two-channel opening was found to be 
an inverse square-root, of much slower decay than the 
RMT prediction. In a later work [23|, this behavior was 
related to the correlations between real and imaginary 
parts of the wavefunction, which in general may depend 
on the underlying classical dynamics. 

On the other hand, real- and phase-space localization 
have been detected in closed systems in correspondence 
of the so-called scars [13, [13|. Within that framework, 
the distribution of the intensities on an unstable peri¬ 
odic orbit was found to decay slower than the RMT- 
expected (10 . due to a phenomenon of constructive in¬ 
terference. This is not our case: in order to rule out 
scarring, we place our probe states away from periodic 
orbits. Still, the localization found for weak coupling to 
the opening does hold in the semiclassical limit, which 
makes us think of a classical effect. 

Successively, we follow the evolution of the wavefunc¬ 
tion statistics of the quantum map for strong coupling to 
the opening. As a result, the intensity distribution be¬ 
comes separated into several long-lived- and a few short¬ 
lived eigenstates. We show that their intensities are pro¬ 
portional to their decay rates, arguing that this quantum 
effect can be explained with the existing theories on res¬ 
onance trapping [23,[13|. particular, the intensities of 
the long-lived states depend on the escape rates through 
a linear relation akin to Breit-Wigner law 0. 

In the second part of the paper we perform analogous 
simulations on the classical cat map and, by looking at 
the statistics of the eigenfunctions of the classical prop¬ 
agator (Perron-Frobenius operator [l^l), f^e de¬ 

viation from the closed system in all similar to the quan¬ 
tum case for weak coupling to the opening. This obser¬ 
vation corroborates the hypothesis of a classical mecha¬ 
nism behind localization, in this regime. On the contrary, 
we show that a strongly-coupled opening does not result 
in resonance-trapping, which makes the classical setting 
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substantially different from the quantum, in this regime. 

The paper is organized as follows: in section III Al we 
calculate the deviation of the wavefunction statistics from 
an exponential distribution due to a single-channel open¬ 
ing by using first-order perturbation theory. In section 
III Bl we verify the theoretical expectation using random 
Hamiltonians drawn from the Gaussian unitary ensem¬ 
ble (GUE) 0, and successively on the eigenfunctions of 
the quantized cat map [2^. Section fll Gl deals with the 
strong-coupling regime: we analyze the proportionality 
between escape rates and intensities, while we account for 
the localization patterns of left and right fastest-decaying 
eigenfunctions in section III PI In section IIIII we intro¬ 
duce the Perron-Frobenius operator of the same test-map 
as a classical propagator, and numerically demonstrate 
an analogous deviation from RMT of its eigenfunction 
statistics for both weak and strong couplings to a small 
opening in the phase space. Summary and conclusions 
are given in section lYl 


II. WAVEFUNCTION INTENSITY 
DISTRIBUTION 

A. Theory 

Suppose Hq is a GUE Hamiltonian. Since its eigen¬ 
functions are complex valued, their intensities x = 
|(a|'!/)o)P at a certain state |a) follow the exponential dis¬ 
tribution Q 


P{x) = e-U (1) 

Now we open the system at |a) [3 

H = HQ-i^\a){al ( 2 ) 

and ask how the distribution of intensities z = |(a|' 0 )p 
is changed with respect to the exponential, when E is 
small enough. By using perturbation theory [sol - l^ . we 
expand the amplitudes (a|'!/') in the first order as 

(alV^n) - {a\'ipl) - iV {a\^l) ^ ^ x ■ (3) 

Zyh/n ^p) 
p^n ^ 

Left and right eigenfunctions are in general distinct for 
the non-hermitian operator but they are just the 
complex conjugate of each other in first-order perturba¬ 
tion regime. We recognize two uncorrelated quantities, 
^ whose real and imaginary parts are Gaussian 

distributed , and rj = ’ following 

(4) 


Pi.z) 



P{z) 



FIG. 1. (a) Rescaled sample distributions of the overall wave- 
function intensities P (2 = | (a|'!/’} 1^) in log scale, obtained di¬ 
agonalizing several realizations of a GUE Hamiltonian, for: 
N = 16384 (dots, 8 realizations), N — 4096 (diamonds, 18 
realizations), N = 200 (squares, 600 realizations), and loss 
parameter F = 0.5; solid and dashed lines are the theoretical 
expectation and the exponential distribution 0 , respec¬ 
tively. (b) The same analysis with the quantum cat map (1141) : 
N = 4096 (diamonds, 28 realizations), N = 200 (squares, 600 
realizations), and F = 1. 


bution of the variable 2 ; = |(a|' 0 )p = namely 

P{z) = J d^drjS {z-e- PoiOPliv) 


^ / _ 

TT J (1-b 7^r7^)^(l-b r^?7^) ’ 


(5) 


where Po{^) oc e . We immediately see that its expec¬ 
tation value 


{z) = 




( 6 ) 


always exceeds unity, meaning the opening produces a 
longer tail, and therefore a certain amount of localization 
of the probability density takes place. 


B. Numerical tests 


with 7 = AEtt ^, and AE ave rag e level spacing of PIq We now verify the theoretical intensity distribution (0) 
(derivation in Appendix El and [30j). We seek the distri- first by diagonalizing multiple realizations of the non- 
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Hermitian Hamiltonian ([2]), where both Hq and the am¬ 
plitudes (a|'0°) are drawn from the Gaussian unitary en¬ 
semble (GUE). The resulting probability distribution for 
the wavefunction intensities | (a|'0) p in first-order per¬ 
turbation regime agrees with the expression ([5]) as shown 
in the example of Fig.llja). The dimension of the Hilbert 
space chosen ranges from N = 200 to iV = 16384, sug¬ 
gesting that the result holds in the semiclassical Imit. We 
will go back to this issue in sec. IHII 

Figure [Ijb) shows that our prediction for a perturbed 
GUE Hamiltonian also fits the distribution of the wave- 
function intensities of the quantized kicked cat map with 
a small opening. The classical evolution of the cat map 
reads Hill 

F,=Foo (7) 

with 

and 

A/, :(«:) = ("--■“P"!’)) modi (9) 

The quantization of the map is given by [1^ [s^ 

Ue = UoV, ( 10 ) 


where 

{qj\Uo\qk) = 

and 

{qj\Ve\qk) = ^QNi{-e cos 27 vp„,+ 2 ^{qi-qk)pm) ^ 

Pm 


( 11 ) 

( 12 ) 


The quantization of the linear map ([HI) is known to 
possess pseudo-symmetries [s^ that make the spectral 
statistics deviate from the Gircular Unitary Ensemble 
(GUE), hence the use of the perturbation dO]) to restore 
the RMT behavior. Here the opening is a minimum- 
uncertainty Gaussian wavepacket 


(lla) 



1/4 

Q-iq-qo)^/2h+ipo{q-qo)/h ^ 


(13) 


whose center (goiPo) is chosen at random on the unit 
torus (the scar at the origin Bill is carefully avoided). 
The non-uni tary pr opagator is realized by replacing 



= (l - § |fl) («l) Uc. (14) 

All the steps of the derivation of Eq. ([5]) would still hold 
in this case, except for Eq. 0, since the quasienergies of 
the cat map follow the statistics of the CUE, instead of 
the GUE’s. Still, both are asymptotically equivalent for 
fV oo Q . In our simulations we alternatively set N = 
200 and 4096, and produce an ensemble statistics of over 
10® states, by repeatedly diagonalizing the matrix m 
over different values of the kick strength e, chosen at 
random within the range [0.1,0.2]. 



FIG. 2. (a) Overlaps between the opening and the left eigen¬ 
functions of the quantized cat map vs. decay rates 7, showing 
the resonance trapping effect, (b) The linear part of the data 
is well described by Eq. (fT51) . here F = 2%/^. 


C. Strong coupling to the opening 


When we further increase the coupling F in the propa¬ 
gator da, the curve ([S]) no longer fits the numerical data, 
as we leave the perturbation regime. A few short-lived 
left eigenstates are localized on the opening, while the 
rest are characterized by intensity suppression together 
with small decay rates. We will clarify the localization 
patterns of left and right eigenfunctions in section III PI 
while we focus for the moment on the left ones. The 
overlaps between the open region and the eigenstates are 
presented as a function of the decay rate in Fig. HJa), 
reminiscent of the so-called ‘resonance trapping’ effect 
[ 2 ^ [13, Issj - lil] , whose main results we summarize as fol¬ 
lows. 

Consider the complex eigenvalues of i7, En — i^n, 7n 
being the decay rates. It has been observed and ex¬ 
plained [ 23 , that when the overall loss ic = 7 „ is 

greater than the energy range AE where the levels are lo¬ 
cated, there exists one particularly short-lived state j'i/'i), 
having decay rate 71 = in — 0{AE/w), while the rest 
of the modes have = 0{AE/w), so that they are 
‘trapped’ near the real axis, although still complex val- 
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ued. We will now use this property together with a. P — Q 
projection formalism to explain the linear dependence 
of the intensities |((^„^i|a)p on the decay rates in this 
regime. Let PHP be the projection of the Hamiltonian 
onto the fast-decaying state, P = |'0i)((/)i|, and QHQ the 
projection on the remaining states, Q = 

We first write an eigenvalue of QHQ as 

{Ej - i-fj) = {(j)j\QHQ\il;j) 

= {(l)j\QHoQ\tjjj) - i^{(l)j\Q\a){a\Q\^jll5) 



On the other hand, we know that QHQ is almost hermi- 
tian, so that, to a very good approximation, 

iQl^P,))^ = {^,\Q. (16) 

We can now recognize the eigenvalues as 

Ej-ijj ~ {(l)j\Ho\(l)j) - i^\{a\(l)j)\^, (17) 

where the first term is the expectation value of a hermi- 
tian operator, hence a real number, and therefore 



7, = ^|(a|<^,)^ (18) 

so that we ‘return’ to a Breit-Wigner kind-of law, as ver¬ 
ified in Fig. [2Ib) for the simulations of the cat map. 


FIG. 3. (Top) Husimi distribution of the fastest decaying 
eigenstate of the quantized cat map, and (bottom) the left 
and right eigenfunctions the same mode, when the map is 
open with loss parameter F = 2'\/2. The opening is placed 
exactly at the position where the Husimi distribution of the 
left eigenfunction is localized. 


D. Left and right eigenfunctions 

We notice that in the strong-coupling regime the left 
and right eigenfunctions of the propagator ca are well 
distinct. In particular, we show in Fig. [3] the Husimi 
distributions of the fastest-decaying eigenstates, whose 
left eigenfunctions only are supported on the opening. 
This is explained as follows. 

The discrete-time evolution operator m is indeed 
split into unitary evolution t/o, and a projection describ¬ 
ing the opening, P^ = 1 — |a)(a|, so that U = PoUq. 
Given the eigenvalue Xj and its eigenfunctions {(j)j\ and 

{(j)j\PoUo\x) = X*{(j)j\x) 

{xlPoUolipj) = \j{x\ipj), (19) 

the projection Po acts first on the left eigenfunction, so 
that in order for the loss to be maximal the amplitudes 
{(f>j\x) should be supported on the opening, in our case 
the coherent state |a). On the other hand, the unitary 
propagator Uq acts first on the right eigenfunction |'0j): 
in one time step we approximate the quantum evolution 
with the classical map F{x), and 

\{x\Uomf^\E{i:,ix))f, ( 20 ) 

so that loss/decay rate are highest if \tpj) is supported 

on the classical preimage of the opening, F~^{o) 1^ 


(Fig. [3]). The localization patterns of left and right eigen¬ 
functions will differ most when they occur where the sys¬ 
tem is more sensitive to initial conditions, typically away 
from fixed points or stable/unstable manifolds of the clas¬ 
sical map. 

In general the outcomes depend on how the propaga¬ 
tion and the loss are arranged, which is usually U = PqUq 
as in our model, but can be inverted sometimes [^ . 


III. CLASSICAL SYSTEM 

In this section we consider a classical chaotic map with 
a small opening, again looking for deviations from RMT 
of the sample distributions of the wavefunction intensi¬ 
ties, properly defined. The idea is to fit the numerical 
data with analytic formulae obtained equivalently to ([5]) 
in perturbation regime, to then extend the analysis to a 
strongly-coupled opening, as done in the quantum set¬ 
ting. 

Using the density operator p, the wavefunction inten¬ 
sities in the quantum regime can be written as 

\{am^ = {a\p\a). (21) 

Here p obeys the Liouville-von Neumann equation 

ihdtp=[H,p\, (22) 
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whose classical analog is 

dtp={H,p}. (23) 

The classical Liouville propagator can be written as 

Uh = (24) 

where L = {H, •} is the Liouville differential operator. In 
the Hamiltonian case U = —L, and therefore the evolu¬ 
tion (|24l) is unitary. The classical evolution operator is 
supported on a space of generalized functions, and its the 
spectrum has a discrete and a continuous part (Stone’s 
theorem); all the eigenfrequencies lie on the unit circle. 
In particular, ergodic and mixing systems only have one 
isolated eigenvalue, e“° = 1, while the rest of the spec¬ 
trum is continuous 

In reality every system experiences noise, coming for 
example from uncertainties or roundoff errors. How¬ 
ever small, noise breaks unitarity and changes the spec¬ 
trum of the Liouville propagator, from continuous to dis¬ 
crete [s^. The (‘leading’) unit eigenvalue is still there, 
but the rest of the spectrum moves inside the unit cir¬ 
cle. In a closed system, the ground-state eigenfunction of 
eigenvalue equal to unity (natural measure) is real and 
positive definite, the density to which all initial condi¬ 
tions asymptotically converge. The other eigenfunctions 
are in general complex and called ‘relaxation modes’, as 
they are associated with the decay of correlations 

(gim) ='^e-'^"{9\Pn){Pn\f) (25) 

The classical-to-quantum correspondence was studied 
by Fishman and coworkers [5l|, [s^, who found that the 
formal solution to the classical Liouville equation, called 
Perron-Frobenius operator [here x = {q,p)] 

{C* o p) (x) = f dxo 6{x - f\xo)) pixo,0), (26) 
JM 

when discretized, effectively behaves like the weakly noisy 
operator, and has the same spectrum as the quantum 
propagator of the Wigner function in the classical limit 
(a similar result was shown in [H). 

Based on that, we can say that the noise introduced 
by the discretization washes out the fine details of the 
chaotic dynamics, and makes the random-wave assump¬ 
tion hold for the eigenfunctions of (l26l) . These are com¬ 
plex valued (in the phase space), and therefore their 
squared magnitudes (‘intensities’) \pn{x)\^ follow a 
distribution. Ideally, the classical limit of the minimum- 
uncertainty wave packet would correspond to just one cell 
of the phase-space discretization. Here we want to repeat 
the analysis carried in the quantum setting, and appre¬ 
ciate the difference in the statistics of the eigenfunctions 
from the closed to the open system. We believe this is 
done most effectively by taking the sum of the square 
magnitudes over a small phase-space interval, as 

in= dxo dx\pnix)\'^S{x - Xq), (27) 

Jmo J 



q q 

FIG. 4. Absolute values of (a) a fast- and (b) a slow- 
decaying eigenfunctions of the evolution operator (I26II for 
the open cat map 0, with the hole located in the square 
[—0.4, —0.3) X [—0.2, —0.1[, obtained diagonalizing a 10"’^ x 10^ 
discretization (I29II . The manifold structure of the cat map is 
shown in (b). 

that is the overlap of \pn{x)\^ with a delta function [clas¬ 
sical limit of the coherent state |a) of Eq. m] sup¬ 
ported on the probing region Mo- The quantum ana¬ 
log of Eq. (123) would be |(a'|'0)|^, over a num¬ 
ber of probe states. In that case, the probability den¬ 
sity P {J2a 'unperturbed system is a 
distribution with M degrees of freedom, 

>M/2-lg-{/2 

Pm{0 - p (iM) 2^/2 ’ 

which becomes a Gaussian as M —)> oo. 

We then perform numerical simulations on the classical 
cat map ©: the Perron-Frobenius operator is discretized 
with Ulam method 

= Ta^ / / dy5{y-f{x)). (29) 

The entries [C]ij are estimated using a straightforward 
Monte Carlo technique , based on counting how many 
trajectories starting from each Mj land in Mi- A partial 
opening is realized by randomly decreasing the number of 
trajectories that start from the hole, which overall covers 
a tiny 1% of the available phase space. The 10^ x 10^ 
matrix is then diagonalized. Fig. SKa) shows a fast- 
decaying eigenfunction peaked in correspondence of the 
hole, as an extreme case of density enhancement at the 
opening. We then measure the statistics of the intensi¬ 
ties (ITTl) in both the closed and open systems (Fig. [5|): 
while the sample taken from the closed system agrees 
with the law (|2^ (M is fitted from the data), the ‘in¬ 
tensities’ on the opening exhibit a longer tail, like in the 
quantum regime. We qualitatively account for this ob¬ 
servation by performing the convolution (0) on the un¬ 
perturbed distribution (E51) . this time in M degrees of 
freedom, 

P{z) oc f r^-MrdpS (z - 

( 30 ) 
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P{z) 



of the eigenstates of the linear operators we are using. As 
a consequence, some of these eigenstates also decay fast 
and are more concentrated on the opening and its preim¬ 
ages For weak coupling, both classical and quantum 
simulations fit this physical picture, and behave likewise. 
Moreover, the calculated deviations of the intensity dis¬ 
tributions from the RMT results all rely on perturbation 
theory, which can be applied to any linear operator with 
a discrete, non-degenerate spectrum. That is the case 
for both the quantum Hamiltonian/propagator and the 
discretized classical evolution operator. 

On the other hand, classical and quantum systems be¬ 
have differently when strongly coupled to the opening, 
the latter only displaying resonance trapping, while the 
former not showing any signatures of mode interaction. 


FIG. 5. Sample distributions of the intensities ^ [given by 
Eq. (I27|l ] of the eigenfunctions of the Perron-Frobenius oper¬ 
ator for the cat map 0: (diamonds) closed system; (filled 
dots) 25% partial opening; (squares) 50% partial opening; 
(empty dots) 75% partial opening; (dashed curve) Eq. (1281) 
rescaled to the data set, where the number of degrees of free¬ 
dom has been fitted from the data to M — 46; (solid curve) 
Eq. (13211 with M = 46, while f = 1 is fitted from the data. 
Inset: total opening at the same location; the peak at the tail 
represents the instantaneous-decay states. 


where |p„(a:)p (here x is discretized by our 

grid), while the perturbation 77 follows [s^ 


Pm(t?) oc 




l + M/2 


(31) 


The outcome is 


Pm{z) = Cm 



(1 -h + y2^2)l-hM/2 


(32) 


where 7 = tt \ Cm = while f is fitted 

from the sample distribution. Fig. [5] also shows two sam¬ 
ple distributions of the intensities obtained for stronger 
couplings to the opening, away from the perturbation 
regime: importantly, the trend of a flatter curve with a 
longer tail stays qualitatively the same, indicating an in¬ 
creasing number of fast-decaying states. A total ope ning 
introduces a number of instantaneous-decay states 
that completely localize on the hole. That generates a 
peak at the very tail of the sample distribution, whose 
shape remains otherwise qualitatively the same as for the 
partial openings (inset of Fig.jSj note the scale). As seen, 
the quantum system in the same regime behaves differ¬ 
ently, as the states that do not decay instantaneously are 
instead long-lived, and the overall intensity distribution 
is consistent with the resonance-trapping picture. 

We may now give an interpretation of our findings. An 
open system, be it classical or quantum, must allow for 
some fast-decaying initial conditions, among the others. 
Densities and wave functions must be expressible in terms 


IV. SUMMARY 

We have shown that: 

i) The overall wavefunction intensity distribution of 
a random (GUE) Hamiltonian and a quantized chaotic 
map deviates from the predictions of RMT, when a 
weakly-coupled, single-channel opening is introduced. 
The result holds in the semiclassical limit. 

ii) By further increasing the coupling to the open chan¬ 
nel in our model, few states localize on the opening par¬ 
ticularly strongly and decay fast, while the rest show the 
opposite behavior: slow decay together with intensity 
suppression at the opening. Using well-known results in 
the context of the resonance trapping effect, we derived 
a linear relation between the intensities of the long-lived 
states and their decay rates, similar to Breit-Wigner law. 
In this framework, we also showed that the difference in 
the localization patterns between fast-decaying left and 
right eigenfunctions can be recognized as an artifact, in¬ 
herent of the construction of open quantum maps. 

in) Analogous simulations of the discretized classical 
evolution operator result in a deviation of the intensity 
distribution from the RMT expectations akin to what is 
observed in the quantum setting, when the coupling to 
the opening is weak enough for perturbation theory to be 
valid. A stronger coupling to the opening increases the 
number of fast-decaying states, so as to obtain a longer- 
tailed intensity distribution, very different from the res¬ 
onance trapping observed in the quantum simulations. 
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Appendix A: Derivation of Eq. (1^ 


We start from the joint probability distribution [s^ of 

Ti — W and d — W 1^1^ 

'/ — l^p^n E^-Ep “ l^p^n (_E„-Ep)2 i 


PivX) oc 


(1 + 7^77^)^ _ 

^2+3M/2 ® 


{l+lW) 
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(Al) 


with M number of degrees of freedom. We simply in¬ 
tegrate over C to obtain the distribution of 77 , in one 
[Eq. dU] or M [Eq. (l3T]) ] degrees of freedom. 
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